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The  p-verslon  of  the  finite  element  method  la  a  new  approach  to  finite  ele¬ 
ment  analysis  which  has  been  demonstrated  to  lead  to  significant  corny- rational 
savings,  often  by  orders  of  magnitude  (this  approach  was  formerly  called  the 
constraint  method;  the  new  term  p-version  is  more  descriptive) .  conventional 
approaches  (called  the  h-verslon)  generally  employ  low  order  polynomials  a? 
basis  functions.  Accuracy  la.ecM.eved  by  snltably  raffelag  the' approximating 
mash.  The  p-version  mass  polynomials  of  arbitrary  or&&  P/Jp2  fox  P*°h-et.a  in 
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plane  elasticity  where  CO  continuity  is  required  and  polynomials  of  order  p  >)' 
for  problems  in  plate  bending  where  Cl  continuity  is  required : 

Hierarchic  elements  which  implement  the  p-version  efficiently  are  used  to¬ 
gether  with  precomputed  arrays  of  elemental  stiffness  matrices. 

Some  important  research  results  are: 

1.  In  polygonal  regions,  under  conditions  usually  satisfied  in  practice, 
if  the  h~version  converges  with  order  of  error  in  energy  norm  0(1/Na)u 
then  the  p-version  converges  with  order  of  error  in  energy  norm  0(1 /n  j  , 
where  N  is  the  number  of  degrees  of  freedom.  This  applies  to  planar 
problems  which  require  either  CO  or  Cl  global  continuity. 

2.  Hierarchic  CO  elements  for  triangles  or  rectangles  have  been  developed 

and  implemented  for  p  >_  2.  Hierarchic  Cl  triangular  elements  have  been 
developed  and  implemented  for  p  5.  Hierarchic  CO  solid  elements  have 
been  developed  for  p  2  for  bricks,  tetrahedra,  triangular  prisms,  and 
rectangular  pyramids. 

3.  Modified  Bernstein  polynomials  have  been  constructed  over  triangles  whict 
provide  a  smooth  approximation  to  functions  in  H^.  These  polynomials  an 
used  to  prove  that  the  p-verslon  of  the  finite  element  method  converges 
in  the  Cl  case  (plate  bending). 

4.  Explicit  mappings  have  been  constructed  which  map  triangles  with  one 
curved  side  Into  the  standard  triangle.  These  mappings  have  the  propert] 
that  the  resulting  elemental  stiffness  matrices  can  be  integrated  in 
closed  form.  The  need  for  numerical  quadrature,  which  may  be  inefficien 
in  the  p-version,  is  thereby  obviated.  The  curved  side  can  be  either 
parabolic  or  elliptic.  The  newly  constructed  mappings  permit  the  p -^ver¬ 
sion  to  be  applied  to  domains  with  curved  boundaries  of  specified  shapes, 
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SUMMARY 


The  conventional  approach  to  finite  element  analysis  is  called  the 
h-version  of  the  finite  element  method .  In  this  approach  the  domain  is  sub¬ 
divided  into  smaller  elements  (usually  triangles  or  rectangles  in  the  plane, 
and  bricks  or  tetrahedra  in  three  dimensions)*  Accuracy  is  obtained  by  fixing 
the  degree  p  of  the  local  polynomial  basis  function  (p  is  usually  taken  to 
be  2  or  3)  and  allowing  h,  the  maximum  diameter  of  the  elements,  to  go  to 
zero*  This  is  in  constrast  to  the  constraint  method  for  finite  elements, 
now  called  the  p-version  of  the  finite  element  method »  in  which  the  subdi¬ 
vision  is  kept  fixed  but  p  is  allowed  to  go  to  infinity.  In  order  to  imple¬ 
ment  the  p-version  efficiently  new  families  of  finite  elements  are  required 
with  the  property  that  as  much  computation  as  possible  can  be  retained  from 
the  p-th  degree  approximation  when  passing  to  the  (p+l)st  degree  approxima¬ 
tion.  Such  families  of  finite  elements  are  called  hierarchic*  When  using 
such  families  the  elemental  stiffness  matrix  corresponding  to  a  pth  degree 
approximation  is  a  submatrix  of  the  elemental  stiffness  matrix  corresponding 
to  a  (p+l)st  degree  approximation. 

The  p-version  uses  polynomials  of  degree  p  2  for  problems  which  require 
C°  continuity  (e.g.  planar  elasticity)  and  polynomials  of  degree  p  5  for 
problems  which  require  continuity  (e.g.  plate  bending). 

Some  significant  research  results  which  were  obtained  are  the  following: 

1.  Rate  of  convergence  of  the  p-version  in  the  case  (singularity 
problem).  Let  ft  be  a  bounded  polygonal  domain  in  the  plane  and  consider  the 
model  problem 


-Au  +  u  *  f  on  ft 


-2- 


where  ru  ■  u  or  ru  ■  3u/3n.  Let  ct^  be  the  angle  made  by  3  ft  at  a  vertex  A^, 

1  i  i  <  Let  Up  be  the  solution  of  (1)  using  the  p-version,  let  u^  be 
the  solution  of  (1)  using  the  h-version,  and  let  uQ  be  the  true  solution  to 
(1)  In  the  weak  sense.  Suppose  N  is  the  number  of  degrees  of  freedom  in  the 
two  finite  element  solutions  Up»  u^.  Then  under  conditions  usually  satisfied 
in  practical  applications,  we  have 


-iy 

1  V"o*  •l.n  m  0(N  2  >  »  I  I^-VqI  li,a  *  0(n"Y) 


where 


Y  *  min  Y^  Y^  «  w/9  ,  1  <_  i  _<  n 


2.  Rate  of  convergence  of  the  p-version  in  the  C'*’  case  (singularity 


problem) . 


Let  £i.  Up,  u^,  o^,  1  i  _<  n,  be  as  before  and  consider  the  model  problem 


A  u  *  f  on  SI 


u  ■  0  du/dn  ■  0  on  3  a  (clamped  edge) 


Let  y^  be  the  smallest  root  of  the  equation 


sin2(Y-l)ai  -  (Y-l)2sin2<»,  -  0 


for  1  <  1  <  n  and  let  y  ■  min  Y. .  Then,  the  conclusion  In  (2)  holds. 

These  two  results  taken  together  show  that  for  two  dimensional  problems 
the  p-version  has  twice  the  rate  of  convergence  of  the  h-version.  This  result 
has  been  demonstrated  computationally  in  extensive  numerical  experimentation. 

3.  Hierarchic  families  of  triangles  and  rectangles  have  been  devel¬ 
oped  and  implemented  for  planar  problems  for  p  2.  Hierarchic  C1  families 


-3- 


*  \ 


of  triangles,  which  use  corrective  rational  functions  to  enforce  conformity 
in  a  minimal  sense,  have  been  developed  and  implemented  for  planar  prob¬ 
lems  for  p  ^  5.  Hierarchic  families  of  solid  elements  for  p  j>  5  have  been 
developed  for  bricks,  tetrahedra,  triangular  prisms,  and  rectangular  pyra¬ 
mids. 

4.  Mappings  of  triangles  with  one  curved  boundary  into  the  standard 
triangle  have  been  constructed  with  the  property  that  the  resulting  ele¬ 
mental  stiffness  matrices  can  be  integrated  in  closed  form.  This  is  impor¬ 
tant  because  in  the  p-version  elements  are  taken  to  be  large  so  that  numeri¬ 
cal  quadrature  may  be  inefficient.  The  mappings  which  have  been  constructed 
are  for  two  cases:  curved  side  parabolic,  straight  sides  arbitrary;  and 
curved  sides  elliptic,  straight  sides  parallel  to  the  axis  of  the  ellipse. 

These  mappings  make  it  possible  for  the  p-version  to  be  applied  when  the 
domain  ft  is  not  necessarily  polygonal. 

3.  Modified  Bernstein  polynomials  have  been  developed.  These  polynomials 
are  defined  over  triangles  in  terms  of  natural  coordinates,  and  they  enforce 
C*  continuity  across  boundaries  of  triangles.  The  modified  Bernstein  poly¬ 
nomials  can  be  used  to  prove  the  convergence  of  the  p-verslon  of  the  finite 
element  method  in  the  case. 
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tory,  Wright -Pat ter son  Air  Force  Base, 

Abstract 

With  one  exception,  all  finite  element  software  systems  have  element 
libraries  in  which  the  approximation  properties  of  elements  are  frozen.  The 
user  controls  only  the  number  and  distribution  of  finite  elements.  The  ex¬ 
ception  is  an  experimental  software  system,  developed  at  Washington  Univer¬ 
sity.  This  system,  called  COMET-X,  employs  conforming  elements  based  on 
complete  polynomials  of  arbitrary  order.  The  elements  are  hierarchic,  i.e. 
the  stiffness  matrix  of  each  element  is  embedded  in  the  stiffness  matrices 
of  all  higher  order  elements  of  the  same  kind.  The  user  controls  not  only 
the  number  and  distribution  of  finite  elements  but  their  approximation  prop¬ 
erties  as  well.  Thus  convergence  can  be  achieved  on  fixed  mesh.  This  pro¬ 
vides  for  very  efficient  and  highly  accurate  approximation  and  a  new  method 
for  computing  stress  intensity  factors  in  linear  elastic  fracture  mechanics. 
The  theoretical  developments  are  outlined,  numerical  examples  are  given  and 
the  concept  of  an  advanced  self-adaptive  finite  element  softward  system  is 
presented. 

2.  "The  Constraint  Method  for  Finite  Element  Stress  Analysis”,  by  I.  N. 
Katz,  presented  at  the  National  Bureau  of  Standards,  Applied  Mathe¬ 
matics  Division  on  October  19,  1977. 

Abstract 


In  conventional  approaches  to  finite  element  stress  analysis  accuracy 
is  obtained  by  fixing  the  degree  p  of  the  approximating  polynomial  and  by 
allowing  the  maximum  diameter  h  of  elements  in  the  triangulation  to  approach 
zero.  An  alternate  approach  is  to  fix  the  triangulation  and  to  increase 
the  degrees  of  approximating  polynomails  in  those  elements  where  more  accuracy 
is  required.  In  order  to  implement  the  second  approach  efficiently  it  is 
necessary  to  have  a  family  of  finite  elements  of  arbitrary  polynomial  degree 
p  with  the  property  that  as  much  information  as  possible  can  be  retained  from 
the  pth  degree  approximation  when  computing  the  (p+l)st  degree  approximation. 
Such  a  HIERARCHIC  family  has  been  formulated  with  p  2  for  problems  in  plane 
stress  analysis  and  with  p  5  for  problems  in  plate  bending.  The  family  is 
described  and  numerical  examples  are  presented  which  illustrates  the  efficiency 
of  the  new  method* 

3.  "The  P-version  of  the  Finite  Element  Method",  by  I.  N.  Katz  and  B.  A. 
Szabo,  presented  at  Air-Force  Flight  Dynamics  Laboratory  Wright- 
Patterson  Air  Force  Base  on  April  23,  1981. 

Abstract 


The  theoretical  basis  of  the  p-version  of  the  finite  element  method  has 
been  established  only  quite  recently.  Nevertheless,  the  p-verslon  is  already 
seen  to  be  the  most  promising  approach  for  implementing  adaptivity  in  practi¬ 
cal  computations.  The  main  theorems  establishing  asymptotic  rates  of  conver¬ 
gence  for  the  p-verslon,  some  aspects  of  the  algorithmic  structure  of  p-version 
computer  codes,  numerical  experience  and  a  posteriori  error  estimation  will  be 
discussed  from  the  mathematical  and  engineering  points  of  view. 


